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We describe an approach for characterizing the process performed by a quantum gate using quan- 
tum process tomography, by first modeling the gate in an extended Hilbert space, which includes 
non-qubit degrees of freedom. To prevent unphysical processes from being predicted, present quan- 
tum process tomography procedures incorporate mathematical constraints, which make no assump- 
tions as to the actual physical nature of the system being described. By contrast, the procedure 
presented here assumes a particular class of physical processes, and enforces physicality by fitting 
the data to this model. This allows quantum process tomography to be performed using a smaller 
experimental data set, and produces parameters with a direct physical interpretation. The approach 
is demonstrated by example of mode-matching in an all-optical controlled-NOT gate. The tech- 
niques described are general and could be applied to other optical circuits or quantum computing 
architectures. 
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Quantum information science promises information 
processing and transmission capabilities far beyond what 
is achievable using classical physics. In particular, quan- 
tum computing has the potential to solve problems which 
are intractable on classical computers. Processing quan- 
tum information requires quantum gates designed to im- 
plement unitary transformations on up to a few qubits. 
In practice experimental quantum gates perform a pro- 
cess which approximates the desired unitary operation 
and usually includes some decoherence. Characteriz- 
ing these quantum processes is critical. This can be 
achieved through quantum process tomography (QPT) 
[H [1 [J which expresses the experimental pro- 
cess in terms of a basis of unitary operations which span 
the space of allowed operations. Experimental noise usu- 
ally results in unphysical process reconstructions and a 
maximum-likelihood correction procedure [(J, Q must be 
used to find the nearest physical process. An alterna- 
tive approach, which we present here, is to construct a 
physical model of the experimental gate by extending the 
Hilbert space to include non-qubit degrees of freedom. 
Experimental data is then fitted to this model to infer 
the parameters describing the system. This step replaces 
the standard maximum-likelihood correction procedure. 
The quantum process in the qubit space is inferred by 
tracing out these additional degrees of freedom from the 
model, whereupon mixing and decoherence effects mani- 
fest themselves. 

A quantum computer is a large interferometer 8] , and 
it is likely to be mode-mismatch of interfering modes 
that ultimately limits performance. We are therefore 
motivated to model imperfections in gate performance 
in terms of mode-mismatch. This requires us to consider 
the non-qubit degrees of freedom of physical qubits. This 
approach is most natural for optical gate implemcnta- 
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tions, whereby the spatio-temporal structure of photonic 
qubits must be considered, however, in principle it could 
be applied to quantum gates in any physical architecture. 

Here we consider linear optics quantum computing 
(LOQC) in particular. We construct a model 
which explicitly allows for the effects of mode-mismatch, 
whereby photon indistinguishability is compromised 
within a circuit, thereby undermining interference effects. 
We show that by fitting the parameters in the model to 
experimental data, the mode-matching characteristics of 
an experimental gate can be obtained. From the model, 
simulated measurement probabilities can be determined 
analytically, which are immune from the effects of exper- 
imental noise. This allows QPT to be performed without 
requiring maximum-likelihood correction to ensure phys- 
icality. This procedure could be applied to other quan- 
tum computing architectures by identifying the physical 
processes of importance and constructing a suitable gate 
model. Importantly, this technique reduces the range of 
input states and measurement bases required to recon- 
struct the process, which may be of considerable advan- 
tage for architectures where the full range of measure- 
ments is not possible. 

Several in-principle demonstrations of element ary 
LOQC gates have recently been performed pTol fllL Ha . 
Il3j . We illustrate our techniques by example of an LOQC 
implementation of the CNOT gate 01 > shown in Figure 
^ The gate employs dual-rail logic whereby a qubit is 
encoded across two spatial modes of a single photon. The 
gate is non-deterministic and post-selected upon detec- 
tion of exactly one photon across the control modes and 
one across the target modes. By considering a gate which 
operates using coincidence detection we significantly sim- 
plify the decoherence effects we need to consider. How- 
ever our a ppr oach can be also used to analyze heralded 
gates |l5l Il6j . Our experimental gate construction is 
identical to that reported in • Experimentally, beam- 
jBsyltbteiBj' arc implemented using waveplates and polariz- 
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ing beam displacers, allowing splitting ratios and phase 
delays to be set with a high degree of accuracy. We there- 
fore restrict ourselves to considering the effects of mode- 
mismatch. 




FIG. 1: Schematic of the CNOT gate using beamsplitters 
with reflectivities r\. c and t denote the control and target 
qubits respectively. Modes labeled ' x ' are discarded and serve 
to balance the amplitudes in the different paths. The gate is 
post-selected upon detection of exactly one photon between 
the c-modes and one between the f-modes. r-boxes represent 
mode-mismatch parameters, described in the text. We adopt 
the phase- asymmetric beamsplitter convention, whereby sign- 
inversion takes place upon reflection from the gray beamsplit- 
ter surfaces. 

We introduce a representation for photons which ex- 
plicitly captures their spatio-temporal structure. Specif- 
ically, we represent photons in terms of their wave- 
function across the photon degrees of freedom (e.g. space, 
time, polarization etc.), 
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where ip(ki, ■ ■ ■ , k n ) is the joint wave-function across the 
photon degrees of freedom k\ . . . k n , and aj^ k is the 
single photon creation operator at the corresponding in- 
finitesimal position, time etc. This is a more generalized 
version of the representation adopted in [ffil | when study- 
ing input distinguishability effects in LOQC. 

The CNOT gate described strictly requires complete 
path indistinguishability at all locations in the circuit 
where photonic interactions take place. When pho- 
ton indistinguishability is compromised non-ideal gate 
operation ensues and the gate no longer implements 
the CNOT logical transformation. The introduction of 
such distinguishability is generically referred to as mode- 
mismatch and is one of the major challenges facing exper- 
imental implementations. Mode-mismatch may arise for 
a number of reasons: imperfect spatial overlap between 
photons; imperfect temporal synchronization; differing 
center frequencies or bandwidths; differing polarization; 
or, any other effect which introduces distinguishing in- 
formation between photons. These problems manifest 



themselves not only for separate photons, but also for 
single photons where self-interference takes place between 
different paths. 

We model mode-mismatch by introducing displace- 
ments (i.e. distinguishability) into the photon wave- 
function at different points in the circuit. Figure ^ shows 
the five locations where displacements are introduced, 
labeled r% . . . t§ . These displacements are sufficient to 
model arbitrary mode-mismatch effects (i.e. the intro- 
duction of additional displacements will be redundant). 
When a photon passes through a r-box it undergoes the 
transformation 
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where r m .„ is the displacement introduced at location m 
into the n th photon degree of freedom. We have assumed 
ij)(k\, . . . , k n ) to be Gaussian for simplicity. Arbitrary 
forms for the wave-function could be chosen, however, 
specifically using Gaussians does not affect the generality 
of the model or detract from its predictive power. This 
is because it is the degree of wave-function overlap which 
is of significance, and, regardless of its form, there will 
always be a set of r's corresponding to a given degree of 
overlap. However, depending on the form and variance 
of the wave-function, the magnitude of the r parameters 
corresponding to a given degree of distinguishability will 
change. This complicates a quantitative interpretation 
of the parameters. Instead they are best interpreted in 
terms of their relative magnitude. 

In addition to mode-mismatch, the model accommo- 
dates for distinguishability which arises during state 
preparation (implicitly incorporated into t%, T2 and T3) 
and measurement (t\ and r 5 ). Parameters T\ and r 5 do 
not affect circuit operation when operating in the com- 
putational basis. This is because state preparation (mea- 
surement) in a non-computational basis is equivalent to 
the introduction of beamsplitters before (after) the cir- 
cuit to generate the required superposition. When oper- 
ating in the computational basis these virtual beamsplit- 
ters are not used and therefore do not form interferome- 
ters, rendering terms t\ and T5 irrelevant. Other optical 
circuits could be modeled in a similar way by identifying 
the locations in the circuit where mode-mismatch could 
occur and constructing an appropriate circuit model. 

It is intuitive that, in the context of gate operation, 
mode-mismatch in any single degree of freedom is com- 
pletely equivalent to mode-mismatch in any other single 
degree of freedom. This is because it is the magnitude 
of photon distinguishability, not the degree of freedom in 
which it is introduced, which results in non-ideal gate op- 
eration. However, it is not obvious that mode-mismatch 
which occurs in multiple degrees of freedom is equiva- 
lent to mode-mismatch in a single degree of freedom. We 
now demonstrate that this is the case, and consequently 
a single degree of freedom in photon distinguishability is 
sufficient to model arbitrary mode-mismatch effects. 

We introduce a geometric representation for the mode- 
mismatch parameter space, as shown in Figure|21 We let 
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each axis of the graph represent a particular photon de- 
gree of freedom, and every point on an axis a vector quan- 
tity, representing the mode-mismatch parameters in that 
degree of freedom. From the point of view of gate opera- 
tion, a point on the fci-axis is completely equivalent to the 
corresponding point on the /c2-axis. If mode-mismatch 
occurs in multiple degrees of freedom, the contribution 
from the different degrees of freedom results in a point 
away from the main axes. However, the choice of axes for 
photon degrees of freedom are completely arbitrary and 
therefore the set of equivalent mode-mismatch parame- 
ters must be rotationally invariant in the space of photon 
degrees of freedom. Thus, any point in mode-mismatch 
parameter space can be rotated onto one of the main 
axes, i.e. into a single degree of freedom. Therefore it is 
sufficient to express photons as a weighted integral over a 
single degree of freedom. Equivalently, five scalar quan- 
tities are sufficient to completely characterize the gate's 
mode-mismatch. This geometric argument generalizes to 
higher dimensions. 



(a) 



(b) 
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FIG. 2: Graphical representation of the mode-mismatch pa- 
rameter space, in two degrees of freedom only (fci and /C2). (a) 
Mode-mismatch in a single degree of freedom, represented by 
the point P. (b) Equivalence of mode-mismatch in any sin- 
gle degree of freedom, represented by the points P and Q. 
(c) Equivalence of mode-mismatch in multiple of degrees of 
freedom, M, and a single degree of freedom, N, through ro- 
tational invariance in the choice of axes. 

The inherent equivalence of the mode-mismatch de- 
grees of freedom unfortunately raises obstacles in the ex- 
perimental interpretation of the parameters. Specifically, 
upon inspection of the parameters it is not possible, in 
principle, to determine in which photon degree(s) of free- 
dom the mode-mismatch is occurring. For example, if 
temporal mismatch occurs this will manifest itself in ex- 
actly the same way as spatial mismatch and there is no 
way from the parameters to infer which is taking place. 

Using the gate model we are able to derive analytic ex- 
pressions for arbitrary coincidence measurement expecta- 
tion values given arbitrary input states. We construct an 
8 by 8 matrix M oxp , of experimentally determined co- 
incidence expectation values, and a corresponding ma- 
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of analytic expressions derived from the 



gate model. The rows of the matrices correspond to the 

input states |00), |01), |10), |++), |+- ), |-+), | }, 

and the columns to the corresponding measurement set- 
tings, where |±) = -^(|0) ± |1}). We define the error 
matrix as 



where the absolute value is performed element-wise. 
From M crror we define the maximum and mean errors 
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We also consider the process fidelity Fp @, |(J 0] , defined 

as 



F P = tr (VXA^XBXA 1 ^ 
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where \a and \b ar e the process matrices of the pro- 
cesses being compared. 

We minimize Ci by optimizing across the mode- 
mismatch parameters T1...T5, which we label T m i n . 
These results are substituted back into the gate model 
to generate an optimized model. The parameters T mul 
are tabulated in Table [I] 



Parameter 


T\ 


T2 


T3 


T4 


r 5 


Magnitude 


-0.30 


0.50 


-0.55 


0.10 


-0.45 



TABLE I: Optimized values of n . . . T5 in units of inverse 
photon bandwidth. 

We apply the estimation procedure to the experimental 
CNOT gate and use ^max and <^mcan to compare the gate 
model to the experimental gate. The results are summa- 
rized in Table |nj It is evident that the optimized gate 
model agrees with the experimental data much better 
than the ideal gate model (i.e. where mode- mismatch 
is ignored). The worst-case error observed is on the 
same order as recent maximum-likelihood QPT recon- 
structions performed on the same experimental gate (ap- 
proximately 2.1%) 0. It should be noted that the model 
always predicts pure states in the extended Hilbert space 
and the r-parameters, which completely characterize the 
gate's operation, have a physical interpretation. When 
expectation values are calculated (i.e. photo-detection is 
applied) these additional degrees of freedom are effec- 
tively traced out, which introduces mixture. 



Gate model 


^raax 




Ideal (i.e. no mode-mismatch) 
Optimized (globally) 
Optimized (independently) 


15.25% 
2.97% 
1.94% 


3.26% 
1.32% 
0.67% 



Ddcl 



(3) 



TABLE II: Maximum and mean error between experimental 
[lH and predicted expectation values, using the ideal CNOT 
gate model, optimized gate model, and optimized gate model 
where parameters are estimated independently for each input 
setting. 

We perform QPT on the optimized gate model to 
construct a process matrix, which does not require 
a maximum-likelihood correction procedure since the 
model is inherently physical and self-consistent. From 
this we calculate the process fidelity with the ideal 
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CNOT process, which yields Fp = 0.88, consistent with 
the result produced through maximum-likelihood estima- 
tion of Fp = 0.87. The process fidelity between the pro- 
cess matrices predicted by the two estimation procedures 
is Fp = 0.95, indicating that the processes predicted by 
the two approaches are highly consistent. Thus, using 
a far smaller data set than is required for full QPT (64 
vs. 256 measurements) we predict essentially the same 
process matrix. This significant reduction comes about 
because we only need sufficient data to determine the 
five r-parameters from which any remaining data neces- 
sary for QPT can be inferred. This could be of partic- 
ular importance in architectures where performing all of 
the measurements required for full QPT is prohibitive. 
Instead, the parameters describing the physical system 
could be determined from a smaller set of more accessi- 
ble measurements, from which all of the data necessary 
for QPT can be inferred. 

Next we repeat the estimation procedure where the 
mode-mismatch parameters are estimated independently 
for each input state used, allowing for input state depen- 
dent effects. This yields a significant improvement over 
the model with static parameters, as shown in Table ITU 
This is indicative that input state dependent effects are a 
contributor to loss of gate fidelity. We infer that imper- 
fections in the waveplates and their alignment change the 
mode-matching conditions in the circuit slightly for each 
different input state. Thus, by expanding the assump- 
tions as to the physical nature of the processes taking 
place, one can significantly improve the accuracy of the 
model. 

Experimentally, the target interferometer is assem- 
bled first and therefore ought to exhibit the least mode- 
mismatch. Thus, we expect the T4 parameter to be small- 
est. This is consistent with our estimation of T4. Follow- 
ing this, the control interferometer is assembled. Simi- 
larly, we expect the t\ parameter to be the next smallest. 
Again, this is consistent with our results. Developing 
a complete understanding of circuit operation based on 
these figures, however, is complicated by the fact that 
the quoted parameters are globally optimized, whereas it 
is evident that input state dependent effects are signifi- 
cant. Thus, a complete interpretation would necessitate 
considering the parameters for the complete range of in- 
put states. Nevertheless, the parameters appear to give a 
consistent indication of the relative mode-match quality 
of different parts of the circuit. 

The parameter T2 only influences gate operation 
through interaction between the control and target 
qubits. In other words, in the absence of the control 
qubit, T2 is invisible. In general, parameters which do 
not affect interference will be invisible. For example, if 



we consider a displacement which affects all gate inputs, 
clearly this is non-interferometric and will not affect 2- 
qubit gate operation. However, if we were to embed the 
gate into say a 3-qubit circuit this parameter could well 
become interferometric. More generally, in order to fully 
characterize a given gate we may postulate that it be 
necessary to embed the gate into a higher dimensional 
circuit in order to observe all physical parameters. This 
is an important realization, since in reality we intend to 
operate gates as a part of larger circuits, not in isolation, 
and as such parameters which are invisible when operat- 
ing in isolation could well become critical when operating 
in the context of a larger circuit [Is| . 

We have presented a model for a CNOT gate in the 
LOQC architecture which explicitly allows for the ef- 
fects of mode-mismatch. By fitting experimental data 
to this model we demonstrated that it is possible to infer 
the mode-matching characteristics of experimental gates, 
leading to significantly improved gate models. In the 
worst case, the error margin in the predictions made by 
the optimized gate models were shown to be similar to 
recent maximum-likelihood QPT studies (calculated us- 
ing a smaller data set). Due to the counting statistics 
of the photon sources we expect gate error to exhibit a 
lower bound of approximately 1.5% [l9|. In the case of 
the optimized gate models the error margins were, on 
average, within this lower bound. 

The ability to infer the physical characteristics of a 
gate represents a powerful diagnostic tool. It also al- 
lows us to perform QPT with a significantly reduced 
experimental data set, since only the parameters stip- 
ulated by the physical model need be determined, from 
which remaining measurement probabilities can be in- 
ferred. Because these probabilities are exact, the neces- 
sity for maximum-likelihood correction is mitigated. This 
approach to performing QPT differs from maximum- 
likelihood techniques in that it makes specific physical 
assumptions regarding the nature of the processes taking 
place and models these by expanding the Hilbert space 
to include non-qubit degrees of freedom. We also showed 
that by expanding the physical assumptions the model 
can be significantly improved. The model presented pro- 
duces estimates which are inherently physical, pure and 
self-consistent, and relies on parameters with a physical 
interpretation. The techniques described are non-specific 
and could be applied to other optical circuits, or quan- 
tum computing architectures by first constructing suit- 
able physical models. 
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